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Abstract 

We study universality in the kinetics of spinodal phase separation in unstable thin liquid films, 
via simulations of the thin film equation. It is shown that, in addition to morphology and free 
energy, the number density of local maxima in the film profile can also be used to identify the early, 
late and intermediate stages of spinodal phase separation. A universal curve between the number 
density of local maxima and rescaled time describes the kinetics of the early stage in d = 2, 3. 
The Lifshitz-Slyozov exponent of —1/3 describes the kinetics of the late stage in d = 2 even in the 
absence of coexisting equilibrium phases. 
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The kinetics of spontaneous phase separation, usually referred to as spinodal decomposi- 
tion, is of interest in far-from-equilibrium systems in diverse areas ranging from materials 
science and biological physics to cosmology and astrophysics [lj. Spinodal phase separation 
is characterized by a combination of early, intermediate and late stages. The early stage 
is marked by amplification/relaxation of initial fluctuations and dominance of the fastest 
amplifying mode which emerges as the spinodal wave. Further amplification and saturation 
of the spinodal wave leads to emergence of new phases and begins the late stage of phase 
separation. Domains of different phases grow as larger domains feed on smaller domains via 
diffusion or advection. The intermediate stage features a mix of formation of new domains 
as well as growth of existing domains. The dynamics of the late stage is well understood, 
and has been successfully explained by Lifshitz-Slyozov (LS) theory [2] and its variants. 
In contrast, the dynamics of the early and intermediate stages remains poorly understood. 
For example, it is not easy to study the time- dependence of the domain size as given by 
correlation functions or structure factors, in the absence of well-defined phases. Experimen- 
tal access to these stages suffers from the added difficulty of detecting low-amplitude and 
multi-modal fluctuations. The present letter investigates the early and intermediate stages 
of spinodal phase separation in the technologically and scientifically important system of 
supported unstable thin liquid films [3(. Our study is simulations of a model thin film equa- 
tion [4]. We discuss the unique universal features of these stages and present a useful marker 
to track their kinetics. The LS growth scenario for the thin film system is also validated 
for the late stage. Interestingly, many different growth laws have been reported for the 
late stage Jp, 6]. The similarity between the thin film equation and the Cahn-Hilliard (GH) 
equation [7| indicates that these unique features are likely to be shared by other spinodally 
phase- separating systems. 

As in all spinodal processes, random fluctuations in the free surface of initially flat, 
supported thin liquid films (< 100 nm) grow and evolve into two distinct phases (viz., a low- 
curvature and thinner flat film phase and a thicker high-curvature droplet phase) whenever 
AG shows a minimum and the spinodal parameter, d 2 AG/dh 2 \h=h < (Fig. [1]). Here, AG 
is the excess intermolecular free energy (per unit area), h is the film thickness, and ho is 
the average thickness 8!, l9[ . Notice that the double-tangent construction for AG in Fig. [1] 
shows that the film phase-separates into phases with h = h m and h = 00. This should be 
contrasted with the usual phase-separation problems such as segregation of binary mixtures 
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which are described by a double- well potential, i.e., there are two possible values for the 



equilibrium composition 



10] . Actually, the droplet phase remains bounded due to overall 



volume conservation and so does the maximum thickness. One way to address this is to 
treat the droplet phase as a defect with changing thickness rather than a true equilibrium 
phase. 

We study the thin film equation which models the spatio-temporal evolution of the film's 
surface in supported thin liquid films. This is derived by considering a thickness-dependent 
excess intermolecular energy AG(h), and simplifying the equations of motion under the 
lubrication approximation 4j. The resulting equation can be written as a CH equation with 
a thickness-dependent effective mobility [7J, M(h) = /i 3 /(3/i) (corresponding to Stokes flow 
with no slip). The total free energy is F s [h] = f[AG(h) + •y(Vh) 2 /2]dx = F e + Fj where, F e 
denotes the net excess free energy and Fi denotes the interfacial free energy. In the above 
expressions, 7 and /i refer to surface tension and viscosity of the liquid film, respectively. 
The corresponding CH equation is 



dh „ 

sF = v ' 



h^fdAG 
— V I — 7V 2 /i 



(1) 



3/i \ dh 

where all gradients are taken in the plane of the substrate. Simulations were done for a 
variety of forms of AG(h) to uncover universal features of the early stage, if they exist. 
Here, a long-range van der Waals attraction due to the substrate, and a comparatively 
short-range van der Waals repulsion provided by a nano-coating on the substrate 111], is 
chosen to illustrate the results. The corresponding AG is shown in Fig. [H and has the form 
AG = -A c /12nh 2 - A s /12n(h + 5) 2 . Here, A s and A c (= RA S ) are the effective Hamaker 
constants for the system, which consists of the fluid bounding the film from the top, film 
fluid and a solid substrate (s) or coating material (c). The thickness of the nano-coating is S. 
Mean film thicknesses falling in the spinodally unstable regime, given by d 2 AG/ <9/i 2 1^=^ < 0, 
are also shown in Fig. [TJ 

Equation ([!]) is solved in the following non-dimensional form to reduce parameters: 



OH 



ff3v , 2nh 2 dAG 



V 2 H 



(2) 



\A S \ dH 

Here, H is the non-dimensional local film thickness scaled with the mean thickness ho] the 
coordinates along the substrate are scaled with the characteristic length-scale for the van der 
Waals case (2ti^ / \ A s \) 1 ^ 2 h^; and non-dimensional time T is scaled with (12n 2 ^h^/ 'A 2 ). 



The excess energy term is completely non-dimensionalized as 



2-nhl dAG 1 f 1 - R R 



(3) 



A s | dH "3 [(H + Df 



where D = 8 /ho is the non-dimensional coating thickness. The linear stability analysis of 



We numerically solve Eq. (J2]) in d = 2, 3 starting with an initial small-amplitude (~ 0.01) 
random perturbation about the mean film thickness H = 1. In d = 2, the system size is uLm 
(n ranges from 16 to several thousands). Periodic boundary conditions are applied at the 
lateral ends. A 64-point grid per Lm was found to be sufficient when central differencing in 
space with half- node interpolation was combined with Gear's algorithm for time-marching, 
which is especially suitable for stiff equations. The parameters D and R were chosen so 
that the film is spinodally unstable at H — 1. An increase in D represents a corresponding 
decrease in the dimensional film thickness (ho) for a fixed coating thickness (8). 

A morphology-based classification of phase separation is shown in Fig. [2j Curved-droplet 
defects and the flat-film phase can be clearly seen in the late stage (middle frame). They are 
absent in the early stage (top frame), where the film surface is characterized by a fluctuating 
wave of increasing wavelength, a result of relaxation of stable modes and growing dominance 
of the spinodal wave. The intermediate stage (bottom frame) shows a mix of developing 
and fully-developed defects. The effect of the spinodal wave can be clearly seen in the 
equi-spaced location and number of defects (~ 10 defects in 10L M ). 

Figure [3] shows the decrease in number of hills or defects (local maxima in film profiles) 
with time. It shows three distinct stages: an early stage with exponent —1/4, a late LS stage 
with exponent —1/3, and an intermediate stage. These stages coincide with those found by 
tracking the morphology and validate the number of hills as a good marker of kinetics. The 
early stage ends when the number becomes comparable to that for the spinodal wave, shown 
as a horizontal line in Fig. |3j The late stage is characterized by the start of —1/3 slope. The 
limits of the intermediate stage can then be defined accordingly. An excellent matching of 
length-scales in the late stages as found by Fourier analysis (solid circles in Fig. [3]) provides 
further support to the choice of number of hills as a marker. As expected, Fourier analysis 
fails to provide any characteristic length-scale in the early stages. Experimental results for 
all stages can be consistently described by pooling early-stage data based on number of hills, 
and late-stage data based on Fourier analysis. Interestingly, linear theory (plotted as dash- 



Eq. (T2]) predicts a dominant spinodal wave of wavelength, Lm = 47r/ 
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dots in Fig. [3]) only explains the beginning of the early stage. So, while the spinodal wave 
itself is given by linear theory, the kinetics leading to its emergence can only be described 
by nonlinear analysis. 

These stages can be readily identified in the evolution of the free energy also (Fig. S]). In 
the early stage, there is not much change in F s , F e and Fi. This stage ends with a sharp 
increase in Fi and a corresponding decrease in F e . In the late stage, there is a smooth 
decrease in both Fi and F e . Fi depends on the sum of squares of local slopes. The random 
perturbations at the start result in a high value of the interfacial energy. The smoothening 
of these leads to a decrease in local slopes and Fi (bottom frame of Fig. HJ). The subsequent 
emergence of the high-curvature defects and their growth lead to increase in F^. On the 
other hand, F e depends on the local thickness and decreases with reduction of thickness for 
films in the spinodal regime. The initial increase in the minimum thickness (inset of Fig. |3J), 
a result of reducing amplitude of stable components of the surface perturbation, leads to 
a slight increase in F e . As the minimum thickness passes through a maximum and starts 
decreasing (inset of Fig. EJ), F e also starts decreasing (bottom frame of Fig. Hj). The total 
free energy decreases all the time, showing the spontaneous nature of the phase separation. 
Again, an analysis of the film profiles shows that the classification based on free energy 
is consistent with the analysis based on the number of hills as well as morphology. We 
remark that Fi is the best marker for defining the intermediate stage. The number of hills 
do not provide such a clear indication of the start and end of the intermediate stage. A 
morphological analysis can demarcate the intermediate stage well but only if it is exhaustive 
and covers the complete profile. 

We have also obtained results for the early stage for several other parameter values (R 
and D), system sizes (uLm) and other force fields, including an uncoated substrate (D = 0) 
and a combination of van der Waals and polar interactions These also show the same 
feature of steady decrease in number of hills (with exponent ~ —1/4, similar to Fig. Ej), 
with results separated only in time scales (results not shown). The effect of system size is 
easily removed by plotting the number density of hills (number of hills per Lm) vs. T rather 
than number of hills. One finds that the kinetics is independent of the system size (results 
not shown). This independence is especially useful for experimentalists who need not worry 
about the lateral size of the sample. 

A master curve to describe all these results will strongly indicate a universal early-stage 
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kinetics. How does one rescale the results to arrive at the master curve? It is known 
that initial fluctuations greatly influence the kinetics for film thicknesses at the edge of the 
spinodal regime or in the defect sensitive spinodal region (DSSR), as opposed to thicknesses in 
the deep inside spinodal region (DISR). 12|. A new evolution coordinate found by dividing 
T by a time, Tj, at which different films are at a given morphological event during their 
evolution, will account for the different initial fluctuations. There can be many such choices 
of events, with the emergence of the flat-film phase being the most obvious one. However, the 
increasing dominance of the repulsive force field slows down the dynamics near this event, 
leading to uncertainty in estimating the corresponding times. Another event, which provides 
a much sharper estimate, is found by tracking the minimum thickness of the evolving film: 
Tj can be chosen as the time at which it reaches its largest value (see inset of Fig. Indeed 
all results, when plotted using this evolution coordinate, collapse on a master curve as shown 
in Fig. [5j This shows that the early stages of phase separation in thin liquid films, as marked 
by the number density of hills, can be represented by a universal curve with exponent —1/4. 
Further confirmation of this universality is provided by results for the d = 3 version of Eq. 
(j2J). They also show the same universality but with an exponent closer to —1/3. This is 
consistent with results from the phase separation of mixtures |1|. The early formation and 
coalescence of domains is known to have a dependence on dimensionality. 

In summary, we conclude that the number density of local maxima (defects) and not the 
usual measures like structure factors and correlation functions can be used to identify the 
early, late and intermediate stages of spinodal phase separation in thin films. The kinetics of 
the early stages for a wide range of potentials and parameters can be described by a universal 
curve between the number density of defects and a rescaled time in d = 2, 3. This process 
has strong similarities with the usual phase-separation processes except for the absence of 
coexisting phases. Still, LS exponent of —1/3 governs the late stages of domain growth. Our 
results in this letter will facilitate the understanding of many experimental results for phase 
separation in thin films. We also hope that our predictions of universality in early-stage 
dynamics will be put to experimental test. 
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Thickness (nm) 



FIG. 1: Variation of the free energy (per unit area) of a film (AG) and force per unit volume 
(d 2 AG/dh 2 ) with film thickness. The parameter values are R = —0.1 and 5 = 5nm. Spinodal 
phase separation occurs for ho > 6.25 nm, where d 2 AG/dh 2 < 0. 
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FIG. 2: Early, late and intermediate stages of phase separation in a supported thin film. The 
dotted, solid and dashed lines refer to non-dimensional times 0.6, 15 and 760 (top), 23000,200000 
and 700000 (middle), 2310,2550 and 2725 (bottom). The parameters are R = -0.1 and D = 0.5. 
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FIG. 3: Variation of the number of hills (defects) with non-dimensional time for a system size 
of 2048 Lm- The filled dots show the corresponding results from Fourier analysis. The dash-dot 
line represents the corresponding linear results. Dotted lines with slopes of —1/4 and —1/3 show 
growth exponents for the early and the late stages, respectively. The parameters are R = —0.1 
and D = 0.5. 
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FIG. 4: Variation of the interfacial free energy (F«), excess free energy (F e ) and total free energy 
(F s ) with non-dimensional time. The bottom frame shows the magnified view in the early stage. 
The parameters are R= —0.1 and D = 0.5. 
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Evolution Coordinate 



FIG. 5: Master curves in d = 2,3 for variation of number density of hills with the evolution 
coordinate. We superpose data for a wide range of potential and parameters, as described in the 
text. The inset shows the variation of minimum thickness with T which is used in defining the 
evolution coordinate. 
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